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ISOMETRIC EMBEDDING VIA STRONGLY SYMMETRIC POSITIVE 

SYSTEMS 

GUI-QIANG CHEN, JEANNE CLELLAND, MARSHALL SLEMROD, DEHUA WANG, 

AND DEANE YANG 


Abstract. We give a new proof for the local existence of a smooth isometric embedding of a 
smooth 3-dimensional Riemannian manifold with nonzero Riemannian curvature tensor into 
6-dimensional Euclidean space. Our proof avoids the sophisticated arguments via microlocal 
analysis used in earlier proofs. 

In Part 1, we prove a local existence theorem for nonlinear, first-order PDE systems that 
do not fall under one of the standard types (elliptic, hyperbolic, and parabolic), under an 
assumption that we call strong symmetric positivity. In Part 2, we apply this theorem to 
prove the local existence result for isometric embeddings. 
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1. Introduction 


Let (M, g ) be an n-dimensional C 00 Riemannian manifold. Recall that a C 00 map y : M —* 
M. N is called an isometric embedding if y is injective and the restriction of the Euclidean metric 
on ~R n to the image y (M) agrees with the metric g on M. In terms of local coordinates 
x = (x 1 ,..., x n ) on M, this is equivalent to the condition that 

(1-1) d i y-d j y = g ij , 1 < i, j < n, 

where g = g %3 dx l dx 3 and <9* denotes -£j. 

In this paper, we study the local isometric embedding problem, which asks whether, given 
a Riemannian manifold (M, g) and a point Xo e M, there exists an isometric embedding of 
some neighborhood of Xo into M N —i.e., whether the PDE system (II. ip has local C 00 solutions 
in some neighborhood of xo. The system (11.11) consists of \n(n + 1) partial differential 
equations (PDEs) for N unknown functions y = (y 1 ,, y N ); thus it is overdetermined when 
N < |n(n +1), underdetermined when N > \n{n + 1), and determined when N = ^n(n + 1). 

The isometric embedding problem has a long and active history. The famous theorem of 
Cartan and Janet (see, e.g., |13| ) guarantees that, when the metric g is real analytic, local 
real analytic solutions to (11.11) always exist in the determined case N = \n{n + 1). In the 
C 00 category, much less is known. Nash |21J proved a global existence theorem in the highly 
underdetermined case N = \n(n + l)(3n + 11). Later, refinements were given by Greene 
[5] and Gunther [6] for the local existence problem that improved the upper bound on the 
embedding dimension to N = |n(n + 1) + n. 

When N = \n{n + 1), known results for g in the C' c0 category are limited to n ^ 4. Most 
research activity has been concentrated on the case n = 2, where local isometric embeddings 
of varying regularity have been shown to exist in a neighborhood of any point Xo 6 M where 
either the Gauss curvature /L(x 0 ) is nonzero, /L(x 0 ) = 0 and V/L(x 0 ) ^ 0, or K (x 0 ) vanishes 
to finite order in certain precise ways (cf. [8] O EH El EEH EE ITS]). For a detailed account, 
see [TO] . 

For n ^ 3, there are fewer results. Bryant, Griffiths, and Yang p] showed that, for n = 3, 
local C 00 isometric embeddings exist in a neighborhood of any point x 0 e M where the 
Einstein tensor has rank greater than 1. Subsequent work was able to relax this restriction 
on the Einstein tensor: In [2D], Nakamura and Maeda extended the existence theorem to a 
neighborhood of any point where the Riemann curvature tensor does not vanish, and in [53], 
Poole extended the existence theorem to a neighborhood of any point where the Riemann 
curvature tensor vanishes but its covariant derivative does not. Meanwhile, for n = 4, the 
results of m. a. and [20] imply that there exists a finite set of algebraic relations among 
the Riemann curvature tensor and its covariant derivatives, with the property that a local 
isometric embedding exists in a neighborhood of any point where these relations do not all 
hold. 

Our main result is a new, simpler proof of the following theorem of Nakamura-Maeda [2D] 
when n = 3 and N = 6 (also see Goodman-Yang [4]): 

Theorem (cf. Theorem (2]). Let ( M,g ) be a C 00 Riemannian manifold of dimension 3; let 
x 0 e M , and suppose that the Riemann curvature tensor R(x 0 ) is nonzero. Then there exists 
a neighborhood O c M o/x 0 for which there is a C 00 isometric embedding y : O —>• K 6 . 

Our proof, like previous ones, uses the Nash-Moser implicit function theorem (cf. Theorem 
lA.ip to obtain a solution. This requires showing that the linearized system has a solution that 
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satisfies certain estimates known as “smooth tame estimates.” (This terminology is due to 
Hamilton; see [7].) The advantage of our approach is that it completely eliminates the need 
for the microlocal analysis and Fourier integral operators used in the proofs of Nakamura- 
Maeda and Goodman-Yang; instead, it is based on Friedrichs’s theory of symmetric positive 
systems. 

Friedrichs [3] introduced the notion of a symmetric positive partial differential operator P 
to study a class of first order linear systems of PDEs 

(1.2) Pv = A l diW + Bv = h 

that do not necessarily fall into one of the standard types (elliptic, hyperbolic, parabolic). He 
proved, under suitable boundary conditions on the domain D, the existence and uniqueness of 
an L 2 (D) solution to the system (11.2)1 . No higher order regularity of solutions is guaranteed, 
even if the functions A 2 , B, and h are G 00 . 

We call a domain that satisfies Friedrichs’s boundary condition P-convex (cf. Definition 
12.3p . Such a domain D has the remarkable property that any solution v to a symmetric 
positive system (II.2p on D is unique in L 2 (D), without assuming any boundary conditions 
on v. This surprising rigidity occurs because a symmetric positive operator P always has a 
subtle type of singularity in the interior of a P-convex domain. In §21 we give a 1-dimensional 
example, where the system reduces to a scalar ODE, that illustrates how this occurs. 

We introduce in this paper a new positivity condition that we call strong symmetric pos¬ 
itivity (cf. Definition 12.ip and prove a local existence and regularity theorem for Erst order 
linear and nonlinear systems satisfying it (cf. Theorem [[]). As the name indicates, this 
condition is a strengthening of Friedrichs’s notion of symmetric positivity. Moser [19] intro¬ 
duced a similar but weaker assumption, closely related to the Legendre-Hadamard condition, 
and proved that any real analytic system of the form (11.21) satisfying this condition on a P- 
convex domain has a unique real analytic solution v. Tso [3D] proved a similar G 00 existence 
theorem on a P-convex domain assuming Moser’s condition, but we believe that his proof ac¬ 
tually requires the stronger assumption of strong symmetric positivity. Both Moser and Tso 
used their results for linear systems to prove analogous perturbation theorems for nonlinear 
strongly symmetric positive systems 

(1.3) <f>(u) = f 

on a domain D c R n , provided that f is sufficiently close to <f>(uo) for a given function uo 
and O is P-convex, where P is the linearization of $ at u (J . 

Our proof of Theorem [2] proceeds in two major steps. In Part 1 (§2HSJ), we establish 
the local solvability of a nonlinear strongly symmetric positive system using the Nash-Moser 
implicit function theorem. In Part 2 (§SHEJ), we show that if the Riemann curvature tensor is 
nonzero at x 0 e M, then there exists an approximate isometric embedding on a neighborhood 
of x 0 where the linearized operator can be made strongly symmetric positive by applying 
a carefully chosen change of variables. This argument consists primarily of linear algebra 
and requires essentially no analysis beyond that required for Part 1. Theorem [2] then follows 
by the smooth tame estimates established in Part 1 and the Nash-Moser implicit function 
theorem. 

The first step requires solving linear strongly symmetric positive systems on a sufficiently 
small, but fixed, neighborhood of a point x 0 in the domain and showing that solutions satisfy 
smooth tame estimates. Surprisingly, Tso’s global existence theorem for strongly symmetric 
positive systems on a P-convex domain does not directly imply a local solvability theorem. 
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This is because there does not necessarily exist a P-convex domain in a neighborhood of a 
given point Xo. This subtle fact is best illustrated by the 1-dimensional example given in 
§3j In we show how this difficulty may be overcome by first restricting the linearized 
system to a sufficiently small neighborhood of x 0 and then extending the restricted system 
to a large ball in M n that is P-convex for the extended system. 

Before proceeding, we recall the following standard notations and facts regarding Sobolev 
spaces on a domain c= M n : 

• The Euclidean norm on vectors or matrices is denoted by | • |, and the £ 30 -norm on 
vectors or matrices is denoted by | ■ |oo- 

• The Sobolev spaces are denoted by 


W k ’ p (Sl) = {ue L p (n) : \\u\\ k>p < 00} 


where || u\ k , P = 2| a |<fc II-D^Hl? is ^ ie Sobolev norm for the multi-index a = (ay, .. 


a n ), and D a u = {dx i r i. u {dx n )an • 

• For p = 2, W k ' 2 (£i) is denoted by H k (Q), with the norm || • ( 1^2 denoted by 

• The C' fc (0)-norm is denoted by 


(fo: 1 ) 1 *! ■■■(dx n ) arL ' 


k 



• The Sobolev embedding theorem [28] implies that H k+m (Q) can be continuously 
embedded into C k (Q) whenever m ^ 1 + [^1; in particular, there exist constants M kl 
depending only on 0 , such that 


(1.4) 


u\\k,co ^ M fc ||w|| fc+ i + [„]. 


Part 1. A Local Existence Theorem for Strongly Symmetric Positive Systems 

2. Strong symmetric positivity 


Let c R n be a bounded, open domain with piecewise smooth boundary and coor¬ 
dinates x = (x 1 ,...,x n ). Let : C x (0,, M s ) —> C' 00 (r2,M s ) be a C 00 , nonlinear first-order 
partial differential operator. Explicitly, for u e C' 00 (f2,M s ), write 


$(u) = F(x,u,Vu), 


where F(x, z, p) = (P 1 (x l , z a ,p °;),..., F s (x l , z a ,p \ l )) is a C 00 , IP-valued function on fl x R s x 
M ns . Given a function f e C c0 (consider the PDE system 


( 2 . 1 ) 


*(u) - f. 


The linearization of $ at the function u 0 e C' a 0 (f2,M ,s ) is the linear first-order partial 
differential operator 3>'(u 0 ) : C' 00 (f 2 ,M ,s ) —* (^(f^lP) defined by 


n 


<&(uo + fv) = y A l diV + Bv 


( 2 . 2 ) 


*' ( u«)v - J t 
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where A 1 , B e C 00 ^, M sxs ) are given by 

V(x) = [(V(x)Q = 

B(x) - [(B(x))“] = 


dF a 

dp\ 

dF a 


(x,u 0 (x), Vuo(x)) 
x,u 0 (x), Vu 0 (x)) 


We will also consider the linear PDE system 

n 

(2.3) 2 A ^ v + 5v = h ’ 

where h e C 00 ^, M s ). 


1=1 


Definition 2.1. The linear partial differential operator (12.2[) is called: 

• symmetric if the matrices A 1 (x),..., A n (x) are symmetric for all x e 12; 

• symmetric positive if it is symmetric and the quadratic form Q o(x) : —* M defined 

by 

n 

(2.4) Qo(x)(f) - y(B(x) + B t (x)-^ a/(x)){ 

i =1 

is positive definite for all x e 12; 

• strongly symmetric positive if it is symmetric positive and the quadratic form Qi(x) : 
M ns —> M defined by 

n 

(2.5) <3i(x)(£i, • • • ,fn) = Yj $ (^*( x ) + %A J (x)) & 

*,3=1 

is positive definite for all x e 12. 

The nonlinear system fl2.ll) is called symmetric (resp., symmetric positive , strongly symmetric 
positive ) at u 0 if the linearization (12.2p of <f> at u 0 is symmetric (resp., symmetric positive, 
strongly symmetric positive). 


Remark 2.2. A few remarks are in order regarding Definition 12.11 

• The quadratic form Q i(x) can be represented by the symmetric ns x ns matrix 


2d i A 1 (x) 

d\ A 2 (x) + d 2 A 1 (x) 


dxA n {x) + d n A\x) 

d\ A 2 (x) + d 2 A l {x) 

2 d 2 A 2 (x) 


d 2 A n (x) + d„A 2 (x) 





diA n (x) + ^(x) 

d 2 A n (x) + d n A 2 [x) 


2 d n A n (x) 


We will use the notation (Qi)y(x) to denote the (i,j )th block of Qi(x): 

(Qi)p(x) = diA°{x) + djA'(x). 


• The positivity of Qi(x) is called the Legendre condition ([T7]. p.10). Moser |19j 
established an existence theorem in the real analytic category assuming the slightly 
weaker Legendre-Hadamard condition (sa. p.ll), which requires only that 

((Qi)^) ab (x)r^v>Aierir/i 2 
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(2.7) 























for all £ e M s , 7] e K n , and some A > 0. But in the C 00 category, the stronger Legendre 
condition is necessary m , bid- 

Definition 2.3. Given a linear strongly symmetric positive first order partial differential 
operator P = A l di + B on a domain hi c R n , the domain 12 is called P-convex if the 
characteristic matrix 

n 

/5( x ) = X^( x ) Ai ( x )> 

i= 1 

where z/(x) = (Vi(x),..., ;y n (x)) denotes the outer unit normal vector to <912 at x e <912, is 
positive def ini te at each point x e <912. 

Tso [30] proved the following: 

Theorem (Theorem 5.1, [30]). Suppose that $(0) = 0 and that the system (12. ip is strongly 
symmetric positive at every (7°° function u in some C 1 -neighborhood of the function uo = 0 
on a domain 12 c M n that is P-convex for the linearization P of Q at uo = 0. Then there 
exist an integer (3 and e > 0 such that, for any f e C' 00 (f2,M s ) with ||f 1^ < e, there exists a 
solution u e C' 00 (f 2 ,K s ) to the nonlinear system (12.ip on 12 . 

Remark 2.4. Note that the condition that a PDE system be symmetric is not an open 
condition with respect to the coefficients. Since the Nash-Moser implicit function theorem 
requires solving the linearized equation not just at uo, but at all u near Uo, it is necessary 
to assume that 4>'(u) is symmetric for all u in some neighborhood of u 0 . The positivity 
conditions, however, are open conditions; hence it suffices to assume that they hold at u 0 . 

Moser [T5] proved this theorem in the case where $ and the function f in equation (12. ip 
are real analytic, under the weaker assumption of symmetric positivity together with the 
Legendre-Hadamard condition (12.71) . Tso [30] stated this theorem assuming these same 
conditions; however, we believe that Tso’s proof, which uses the Garding inequality for non- 
compactly-supported vector-valued functions on the domain hi, is correct only if the stronger 
Legendre condition holds. (See [27] and the discussion at [31].) 

3. A LOCAL EXISTENCE THEOREM FOR STRONGLY SYMMETRIC POSITIVE SYSTEMS 

The goal of Part 1 of this paper is to prove the following local version of Tso’s theorem: 

Theorem 1 . Suppose that the linearization <f>'(u) of <$> is symmetric for all u in some C 1 - 
neighborhood of uo e C' co (f 2 , R s ), and that $'(uo) is strongly symmetric positive at some point 
x 0 e 12. Then there exist a neighborhood 12 0 c 12 o/x 0 , an integer j3, and e > 0 such that, 
for any f e <77°°(L2 0 , M s ) with ||$(u 0 ) — f ||/3 < e, there exists a solution u e C' 00 (r2o ; ^ s ) t° the 
nonlinear system (12. ip on h2 0 - 

We wish to emphasize that Tso’s theorem does not immediately imply the local existence 
result, because strong symmetric positivity on a domain does not necessarily guarantee 
the existence of a P-convex neighborhood of xo- In fact, as we show in the example below, 
in general no such neighborhood exists. 

Example 3.1. Consider the following ODE: 

(3.1) (x — x 0 )u' + bu = h(x). 
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It is straightforward to verify that (13. ip is strongly symmetric positive if b > 1, and an 
interval D = (xi,X 2 ) is P-convex if and only if Xo e (xi,x 2 ), i.e., if and only if the regular 
singular point of this ODE lies in the domain. 

Meanwhile, the general solution of (13.11) is 



which is continuous at x = Xq if and only if C = 0. Thus we see that: 

• The P-convexity condition forces the uniqueness of a C® solution of (13.If) on 0, 
without specifying any initial or boundary data for u. 

• If 0 is not P-convex—i.e., if x 0 ^ D, then the ODE (13. ip has infinitely many solutions 
on O. In this case, P-convexity—and hence uniqueness of the solution—can be 
achieved by extending the domain to one that contains the singular point xq. 

In higher dimensions, a similar phenomenon occurs: Consider the strongly symmetric 
positive linear PDE system (12. 3 p on a domain O c K n , and let xi,X 2 e <9D be located on 
opposite sides of dD, with v = z/(xi) = — z/(x 2 ). In order to have ,d(xi),/3(x 2 ) > 0, the 
matrix zyA*(x) must be positive definite at xi and negative definite at x 2 . Therefore, P- 
convexity requires that each of its eigenvalues must change sign somewhere in the interior of 
D. For n ^ 2, this does not necessarily imply that the system (12.31) has any singular points 
in D, but it is still true that any C 00 solution on D is unique. Moser discusses this in [19]. 
concluding that, “The reason for this strange phenomenon is that usually the conditions [of 
the theorem] imply the presence of a singularity and a solution which remains smooth at the 
singularity is unique.” 

Our proof of Theorem Q] will proceed as follows: without loss of generality, assume that 
u 0 = 0 and $( 0 ) = 0 . 

• In 1 14.11 we restrict the nonlinear system (12.11) to a neighborhood Do c D of x o on 
which the quadratic forms (Q u )o( x ) and (< 3 u)i( x ) associated to any sufficiently small 
function u on D 0 remain sufficiently close to Q 0 (0) and Qi(0), respectively. 

• In 1 14.21 we extend the linear PDE system (12.31) from the domain D 0 to a strongly 
symmetric positive system on all of R ri , where the coefficients satisfy C 1 bounds that 
will be needed later. 

• In 1 14.31 we show that, for sufficiently large R > 0, the ball Br of radius R is P-convex 
for the extended linear system. 

• In 1 14.41 we use the extended linear system on Br to prove the smooth tame estimates 
required to implement a Nash-Moser iteration scheme to solve the nonlinear system 


02.ip on D 0 


Appendix [A] contains the precise statements of the Stein extension theorem [29] and the 
Nash-Moser implicit function theorem [26j that will be used in the proof of Theorem [0 


4. Proof of Theorem Q] 

4.1. Restriction of the nonlinear system to an appropriate neighborhood of Xo. 

Without loss of generality, assume that u 0 = 0, <I>(0) = 0, and x 0 = 0. First, we show how 
to choose an appropriate neighborhood D 0 on which to construct a solution for the system 


m- 
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For ease of notation, set 


B = B( 0), A i = Ah(0), A) = djA\ 0). 

Using Taylor’s theorem with remainder, we can write 

n 

(4.1) B(x) = B + B(x), ^(x) =A i + '£ l x i A)+A i (x), 

3 = 1 

where 5, A* e C' 00 (r2, M sxs ) are such that B vanishes to order 1 and A 1 vanishes to order 2 at 
x = 0. The strong symmetric positivity hypothesis at x = 0 is equivalent to the assumption 
that the quadratic forms Q 0 : M s —>► M, Qi : —> M defined by 

n n 

( 4 . 2 ) Q„K) = y(B + B T -y;y)c, <5i(s,...,eo= 2 ej(y + y)& 

i=l i,j =1 

are positive definite. 

Lemma 4.1. Suppose that <h satisfies the hypotheses of Theorem, Q] at x = 0. Let A 0 , Ai > 0 
denote the minimum eigenvalues of Q 0 and Q 1; respectively, and let B r c: M n denote the ball 
of radius r about x = 0. Then, given real numbers M 0 ,Mi > 1 and 5 > 0, there exist real 
numbers r,p> 0 and an integer a > 0 such that B r cz and, for any u e C' 00 (S r ,M s ) with 
||u|| Q < p, the matrix-valued functions B U ,A Z U e C co (B r ,M. sxs ) associated to the linearization 
of <3> at u may be written as 

n 

5 u (x) = B u + B u (x), A’ u (x) =Al + ^ x j (A u ) i j + i’ u (x), 


where 
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(4.3) 

B\i hpoo ^ 2 ’ 
ir ii 5 

\{Au))-A) U <2, 

IIA II 5 

11 u||l ’ 00< 2M/ 

K 

- A% < 8, 

5 


1! -Bull g® < 2Mo ’ 

Hull 

°’°° < Mo' 


For convenience, we will refer to any function u 6 C co (B r ,'R s ) with ||u|| Q < p as “admissi¬ 
ble.” 


Proof. Choose r > 0 so that the restrictions of B and A 1 to the ball B r of radius r satisfy 


(4.4) 


\B\ 


0,oo 


< 


4M n 


II ^ 4 * || 1,00 < 


4Mp 


I v4.' || o,oo < - 


2 M n 


Then the Sobolev embedding estimate (II.4ft and the smallness of the Taylor remainder terms 
for small p imply that we may choose p and a so that equations (I4.3jl hold. Indeed, we may 
choose any a ^ 3 + [|] and then choose p > 0 accordingly. □ 


In 1 14.21 we will show how to choose the constants 8, M 0 , and M 1 so that the restriction 
of the system (12.11) to the domain O 0 = B r has the property that its linearization at any 
admissible u e C co {B r ,M s ) may be extended to a strongly symmetric positive system on all 
of M n . 












4.2. Extension of the linearized system to R n . We will use Stein’s extension operator 
(cf. Theorem IA.2|) to extend the coefficient matrices in the linearized system (12. 3 p from B r 
to all of R n . First we need the following lemma, which states that the bounding constants 
in this construction are independent of r: 


Lemma 4.2. There exist constants Mk :P , 1 < p < go, 0 < k < go, and extension operators 
£ r . : L l (B r ) —» L 1 (M n ) for all r > 0 such that, for all f e W k,p (B r ), 

||Sr/||fc,p ^ ^fc,p||/||fc,p- 

Proof. Theorem IA.2I guarantees the existence of such constants and an extension operator 
for r = 1; then a straightforward rescaling of the operator and a standard rescaling argument 
shows that these constants are independent of r. □ 


Now, set 

Mo = Mo i0 o, Mi = Mi j00 , 

where M 0jO0 and M l oo are as in Lemma [4.21 Choose 5 > 0 such that, for any matrices Ay 
and B' with 

\B'-B\ X <5, \Ay~A)U<5, 1 

the quadratic forms Q' 0 : R s —> R and Q\ : R” s —> R defined by 


ei© = « T (B' + (B') T -XA'k oiKi, ■■•,«»)= Y, $ W+ A ‘P 


i =1 


*p=l 


are positive definite with minimum eigenvalues greater than or equal to ^A 0 and |Ai, re¬ 
spectively. Then take r > 0 as given by Lemma [4.11 and set O 0 = B r . Henceforth, we will 
restrict the systems (12. ip and (12. 3p and all relevant quantities to B r . 

Next, we construct an extension of the linearized system (12.3p on B r to all of R" in such 
a way that the coefficients of the extended system are bounded in VF fc,:P (R ri ) with respect 
to the W k,p (B r ) norms of the coefficients of the original system on B r . After replacing the 
functions B , A *, and h by their restrictions to B r , define C 00 functions B,A l , and h on R ra 
by 


B(x) = B+ (£ r S)(x), 

n 

(4.5) i f (x) =A i + £ i x j A i j + (£ r i i )(x), 

j= 1 

h(x) = (£ r h)(x). 


Similarly, for any admissible u e C co (B r , R s ), let and B u denote the analogous extensions 
of the functions A l u and B u corresponding to the linearization of <f> at u. Then we have the 
extended linear systems 


(4.6) 2 A l u di\ + £> u v = h 

i =1 


on R re . 
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Proposition 4.3. For any admissible u e C co (B r , R s )? the extended system (14.61) is strongly 
symmetric positive on K". Moreover, for any x e R n , the associated quadratic forms 
(Qu)o( x ) : R s —> R and (Qu)i(x) : R ns —> R defined by 

n 

(Qu)o(x)(0 = ^ T (5 u (x) + P([(x) - 2^u( x ))f. 

(4 -7) _ n _ ^ 

(Qu)l(x)(6, • • ■ ,£n) = 2 (^-4’ u (x) + ^A J u (x)) & 

*,i=i 

have minimum eigenvalues greater than or equal to |A 0 and 4 Ai, respectively. 

Proof. By construction, the functions £,P U and £ r A„ satisfy 

r r 

(4.8) ||£r-B U ||o,oo < 2 > ||£r^4ul|l,00 < II||o,oo < & 

The first and second inequalities in (14.8[) imply that, for all x e R n , we have 

B U (x) P|oo^ \B U 51oo + |£r^u(x)|oo < d, 

|^(x) - AJU < |(A U )} - A}la, + |d i (£ r i ? u )(x)| 00 < 6 , 

and the result follows immediately. □ 


4.3. Boundary conditions on Br for large R. Next, we show that, for R sufficiently 
large, Br is P-convex for the extended linear system (14.60 . 

Proposition 4.4. Let R > 0. For x e 8Br, let z/(x) = (z/i(x),..., z/ n (x)) denote the 
outward-pointing unit normal vector to 8Br at x. Then, for R sufficiently large, the char¬ 
acteristic matrix 

n 

(4.9) 0(x) =2^( x )A l u (x) 

i =1 

is positive definite for all admissible u e C' G0 (P r ,M s ) and x e 8Br. 


Proof. The normal vector to the sphere 8Br is given by 

, , 1 


nv = ^x. 


Thus, we have 


(4.10) 


/J( x ) = s X>vi; 

i =1 


= + 2 rxHM +2^(£^)(x: 

i= 1 i,j=l 4=1 

The first and third terms in equation (14.101) are bounded: 

1 n 

Yx l A l u < nmax{|Ai|oo,..., |A£|oo} , 

• ■ rr\ 


(4.11) 


R 


i =1 




2^(£r4)(x) 


4 = 1 


< n<5, 


to 


OO 














where the second equation in (14.111) follows from the third inequality in (14. 8[) . Meanwhile, 
we claim that the second term in equation (14.10)1 has minimum eigenvalue greater than or 
equal to \R\\. This can be seen as follows: Consider the corresponding quadratic form 
(Qu)I(x) : R s —> R given by 


(Ou)I(x)K) = ^ S « 

*,1=1 


J (:r i x^A n YM. 


Then, by Proposition 14.31 we have 


lWu)!(x)«)i = T5 i(d„) 1 (x)( a: i f,..., a ,”f)i 


2 R 

> Ta, (lx'41 2 + ... + |x-{|») 

= 2jA i(( i i)2 + ... + ( J ,” )2 ) | e |2 


- 4i Alfl2 '^ 2 
= jHA^I 2 . 


Therefore, the minimum eigenvalue of (Qu)I(x) is greater than or equal to Together 

with the inequalities in (14.lip , this implies that, for R sufficiently large, /?(x) is positive 
definite for all x e dB R . □ 


4.4. Application of the Nash-Moser iteration scheme. The final step in the proof of 
Theorem [T] is to apply the Nash-Moser implicit function theorem (cf. Theorem lA.ip . 

Notation 4.5. We will adopt the following conventions: 

• Functions without tildes are taken to be defined on B r , and ||v||* will denote the 
H k -noTm of v e H k (B r ). 

• Functions with tildes are taken to be defined on Br, and || v||fc will denote the H k - norm 
of v e H k (B R ). 

Let E k = H k+1 (B r ,R s ) and F k = H k (B r ,R s ), with the usual iL fc -norms; then we have 
Eao = Eos = C°°(B r ,M. s ). Let D 0 c E 0 denote the ball of radius p > 0 centered at u 0 . 

Smoothing operators S(t) : Eq —*■ E x may be constructed as follows (see., e.g., [T| or C5D- 
First, choose a compactly supported function x e C®(M n ) with x 5 s 0 and y(x) dx = 1. 
For t > 0, define 

Xt(x) = t n x(tx), 

and define S t : L 2 (M n ,M s ) -► C7°°(R n ,R a ) by 

(4u)(x) = f Xt( x — y)u(y) dy. 

JR n 

Then, define S t : E 0 —> E <*, by composing S t with the Stein extension operator £ r : L l (B r ) —> 
L\R n ): For ue£ 0 = H\B r ,R s ), deffiie 

(S t u) = (. S t Z r u)\ Br . 

It is straightforward to show that the operators S t satisfy the required inequalities; see [25j 
for details. 


li 












The fact that <f> is C 2 follows from the fact that F is C 00 , and the bounds (1A.2H follow from 
the Gagliardo-Nirenberg and Sobolev inequalities (see, e.g., 0). To complete the proof, it 
suffices to show that there exists an integer a ^ 0 such that, for any integer m ^ a+ 1, given 
any u e D m , the extended linear system (14.611 on B R corresponding to the linearization of 
(12. ip at u has a unique solution v e H m ~ a (B R ) for any h e H a (B R ), and that the restriction 
v = v| b t satisfies the smooth tame estimates (1A.3H . 

First, because the extended system (14.611 corresponding to a given admissible u e H m (B r ) 
is symmetric positive with coefficient matrices (omitting the subscript u to avoid notational 
clutter) A,...,A n ,B e H mL (B R ) and B R is P-convex for (14.6p . Friedrichs’s theory of 
symmetric positive systems [3] guarantees the existence of a unique solution v e L 2 (B R ). 
Moreover, we can obtain an explicit L 2 bound for v, and hence for v, as follows. Multiply 
the matrix equation (14.6p by v T to obtain the scalar equation 

n 

(4.12) ^ vLT diW + v T P v = v T h. 

i=1 

Then, because A 1 is symmetric and 

n n 

2 d t (v T A*v) = 2 (2v T A i (div) + v T (djA*)v) , 

i=1 i =1 

we can write equation (14.12j) as 

n n 

(4.13) - \ v J (di A)v + v T P v = v T h - | d t (v T Abv) . 

2=1 2=1 

Multiply by 2 and use the fact that v T B v = |v T (F + B J ) v to obtain 

n n 

(4.14) v T (P + B J - Yj = 2v T h - d t (v T A*v) . 

2=1 2=1 

By Proposition 14.31 it follows that 

n 

±A 0 |v| 2 < 2v T h ~Yj d i (v T A*v) 

2=1 

^ ^|v| 2 + — |h| 2 -Ydi( v T A i v > ) . 

4 A o ti ' ' 

Integrate over B R , apply Stokes’ theorem, and use the fact that /3 is positive definite on dB R 
to obtain 

(4.15) ||v[|o ^ C' 0 (A 0 ) 2 ||h||o - y- f (v T /3v) dS < C , 0 (A 0 ) 2 ||h|| 2 , 

J8B r 

where C'(Ao) > 0 is a universal constant depending on Ao- Therefore, the restriction v of v 
to B r satisfies 

(4-16) ||v||o ^ ||v||o < C'o(Ao)||h|| 0 < C , o(Ao)-^o, 2 ||h||o, 

where M 0 2 is as in Lemma [4.21 
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Bounds on the derivatives of v may be computed similarly by differentiation. First, dif¬ 
ferentiate the system (14.61) with respect to ad to obtain 

n 

(4.17) ^ (A 1 <9 2 v + (djA l )di\^ + B djv + ( djB)w = djh. 

7 =1 

Multiply the matrix equation (14.171) by dj~v T to obtain 

n 

( 4 - 18 ) 2 ((^v T )i J d 2 v + (d j v T )(d j A t )d i v'j + (d^)B djw + (djV T )(djB)v = (d i v T )d i h. 
By an argument similar to that above, we can write equation (I4.18P as 

n n 

(4.19) <7,-v T (b + B J — ^ diA^dj-v + 2 s P i (dj\ J )(djA l )d i v 

7=1 7=1 

n 

= 2djV J (djh - (djB)v'j (djV T A'djv) . 

7=1 

Now sum equation (14.19j) from j = 1 to n, and note that the second term can be written as 

n n 

2 2 j (djV^djA^div = (djV T ){djA l + diA 3 ){diw). 

i,j= 1 1,7=1 

Thus the summed equation can be written as 

n n n 

(4.20) 2] d,v T (b + S T - 2] dj-v + ^ (d i v T )(d j M* + diA 3 ){d{v) 

j = 1 7=1 7,j = l 

n / n 

= 2] - (^)v) - 2 & (^v T i l 7jV 

7=1 V *=i 



or, in other words, 

(4.21) 

n n 


2] Qo(djv) + <9i(5iv,..., d n v) = SdjV 1 " 
7=1 7=1 V 



(djB)v^j ~Yjdi (dj 

7 = 1 



By Proposition 14.31 it follows that 


(A 0 + Ai) l^-vl 2 s7 2 ( 25 7V T feh - (fyB)*) ~Yi d i (djV T A i d j v 


7=1 


7=i 


7 = 1 




n A n / 

l(Ao + Ai) 2] l^v| 2 + , 7 2 (l 5 7h| 2 + \SjB\l |v| 2 

i=i ^ A ° + Al ^ i=i V 

n 

~Yi di (dj^A'djv) . 


7=1 
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Integrate over Br, apply Stokes’ theorem, and use the fact that /3 is positive definite on 8Br 
again to obtain 


(4.22) 


l|v||? < Ci( A o, Ai ) 2 (J|hj| 2 + |jv|| 2 || J B|| 2 0O ) - 
<C' 1 (A 0 ,A 1 ) 2 f||h|| 2 + |v||ol|5||? 


I1 p 

w) SXj 


(Aq + Ai; — j 8Br 


(djV T f3dj\) dS 


where Ci(A 0 , Ai) > 0 is a universal constant depending on A 0 and Ai. 
By the Sobolev embedding estimate (11.41) . we have 

II-B||i,oo ^ K||B|| 2+ [|] 

for some constant K\ thus we can write the inequality (j4.22[) as 


v 




Cl 



) 


and hence 


v||i < Ci 


h||i + ||v||o||B 


2 +[ 



Therefore, the restriction v of v to B r satishes 


(4.23) 


v||i ^ ||v||i ^ Ci (J|h||i + ||v||o||B|| 2+ [|]J 

< ^||h||i + |jv||o||B|| 2+ [ f] ) 

^ C'l (j|h||i + ||h||o||B|| 2+ [^]^ 

^ C'l (j|h||i + ||h|| 0 ||u|| 3+ [|]^ , 


where the last inequality follows from the fact that B is a (7°° function of u and its first 
derivatives. 

Successive differentiations of the system (14.6j) produce similar results. To obtain an esti¬ 
mate for ||v||fc, differentiate the system (14.61) k times, with respect to x n ,..., x jk . This yields 
an equation of the form 


(4.24) 2 A< V + 2 + B 4 .. - 


i—1 


q=l 


k —1 / n 

- (4,.*S)v T h + D*-™B I (JTv), 

m—1 \i =1 


-) 


where, on the right-hand side, D m indicates an appropriate differential operator of order m. 
Multiply the matrix equation (14.24)1 by 28^ Jfc v T , rewrite the hrst term and rearrange as in 
the previous cases, so that the left-hand side of equation (14.241) becomes 


(4.25) 4..,,v t (B + B t - 2 <^)<4..„v + 2(4..,, v T ) 2 

i =1 i=l q =1 
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Now sum over j\.... , jk , and note that the second term in (I4.25P can be rearranged as follows 
by using the commutativity of mixed partial derivatives and relabeling as appropriate: 

n k 

Ml.-Uk-l 9=1 

n k 

= 2 V V(d fc . , . v T )(d, Ayb. , . v 

AJ 3q3l—3q—3k Jq ’ m—3q—3k 

= l 9=1 

n k 

= V Y(d k . , . v T )(d ? - a A i + d i A j i)d k . , .V 

/—l /—l ■ 3q3l—3q—3k M ' m—3q—3k 

*Jli-dfe=l 9=1 


= * 2 (4..** T )(4^ + 

—iik = l 

Thus the summed equation can be written as 


3k -1 


V. 


Jl>-jJk = l 


jl,-Jfc = l 


■jkY + k 

2 wL.,3 k -i,^-- 

v) 

‘ ’ 3i,-,3k-i,n > 

31 

vifc—1 = 1 


2d k ■ v T 

3l-3k 

(4..,A - (4..*s)v) 


k- 1 

/ n 


-2 

r k u : v 1 2 D k+1 + D k ~ m B (D m v) 

m—1 

\2 = 1 

n 

\ 


-2 ^ ( 4 .. 
2 = 1 



By Proposition 14.31 the left-hand side of equation (I4.26P is bounded below by 

n 

|(A 0 + k\i) 2 \ d h,..,3<y\ 2 ■ 

jl,—,3k=l 

Tlius, after performing operations similar to those above, we obtain 
(4.27) 


fc-i 


v lll ^ C k (\o, Ai) 2 J ||h [| 2 + ||vQMIIl,oo + y 


m—1 


2 II ^111 + 1-771,00 + IIMll-m,oo 

\j =1 


fe -1 

2 , V II ~ II 2 


^ C'fe(Ao)Ai ) 2 ( ||h||fc + y ||v 

m =0 


Xj 1^111+1-771,00 + 1-^111-771,00 

vi=i 


By the Sobolev embedding estimate H 1.4 h . we have 

|| dP || fc+1—m,oo ^ -7^ A; 11 ^4"^ ||/ c _)_2—m+[^]> ||-S|| fc— 771,00 ^ ^ k || B || fc+ 1 —m+[^] 

for some constant Ah,; thus we can write the inequality (I4.27P as 









and hence 


fc-i 


(4.28) 


11% < Ck | ||h||fc + E w m £ ||^||fc+2-m+[f] + \\B\\ fc+1—m+[y] 


m =0 


Vi=1 


By the Gagliardo-Nirenberg interpolation inequality [2] and the Cauchy-Schwarz inequal¬ 
ity, for 0 < m < k — 1, we have 

|| v||rn ||^+ 2 —m+[ 2 .] < C m (|| v|| 0 \\A J || k+ 2 +[^ ] + ||v||fc-l ||A J ||3+[f ]) , 


|| v|| m ||.B||fc+i_ m +["] ^ C' m (||v||o||i?||fe + l + [B] + ||v||fe_i||5|| 2+ [^]) 
for some constant C m . Substituting into equation (14.2811 . we obtain 


(4.29) 


v 11 1 . ^ Cl ||h|L + llvll 


Y ||^||fc+2+[f] + H-S||fc+l+[f] 


\j =1 


+ V \\k-l 


Y l^lla+lfl + l^b+lf] 


\j = 1 


Now let a ^ 4 + [|]. It follows from the fact that A and B are C 03 functions of u and 
its first derivatives that there exist constants K p and K kp such that, for any u e D a , the 
extended linear system (14.611 corresponding to the linearization of (12.ip at u satisfies 

||5||2+[f] ^ K p , ||A'|| fe+2 +[§], \\B\\ k+1+m ^ K k , p { 1 + |u|| fc+3+ [s]). 

Thus (14.2911 becomes 

(4.30) ||v|| fc < C'l (||h|| fe + lvlollulfc+3+f"] + ||v||fc_i) . 

It then follows by induction (with the inequality (14.23ft as the base case) that 
(4-31) ||v||fc ^ C% (j|% + lvlollulfc+3+fn]) . 

Therefore, the restriction v of v to B r satisfies 


(4.32) 


v||fc < ||v||fc ^ C" (j|h|| fc + lvHollullfc+3+fn]^ 

^ CkM k '2 (j|h| fc + ||v||o||u||fc + 3 +[ ~]) 

^ Ck f||h||fc + |h|| 0 ||u||fc +3+ [«]^ . 


All the hypotheses of Theorem I A. 1 1 have now been verified; thus the conclusion of Theorem 
IA.1I gives the desired solution u e C G0 (5 r , M s ) to the nonlinear system (12.111 on B r . This 
completes the proof of Theorem |T| 


Part 2. Application to Isometric Embedding 

5. Local existence theorems for isometric embedding 

The remainder of this paper will be devoted to giving a new proof, based on Theorem [U 
for the following local existence theorem: 
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Theorem 2. Let (. M,g) be a C 00 Riemannian manifold of dimension n = 2 or n = 3; let 
N = |n(n +1); let x 0 e M, and suppose that the Riemann curvature tensor i?(x 0 ) is nonzero. 
Then there exists a neighborhood Q c M of xo for which there is a C 00 isometric embedding 
y:Q^R N . 

Here we briefly describe our strategy for proving Theorem [2] Let n = 2 or n = 3, and 
let N = \n(n + 1). For convenience, choose local coordinates x = (x 1 ,..., x n ) based at 
xo, so that without loss of generality we may assume that xo = 0. Given a (7°° metric g 
on a neighborhood fl of x = 0, choose a real analytic metric g on 12 that agrees with g to 
sufficiently high order at x = 0. By the Cartan-Janet theorem, there exists a real analytic 
isometric embedding (possibly on a smaller neighborhood) yo : 12 —> R N of (12, g) into R N . 

The linearization of the isometric embedding system (11.11) at yo is a first-order PDE 
system of N equations for the unknown function v : 12 —» R N . This system decomposes into 
a system of n first-order PDEs for the tangential components of v, together with (N — n) 
equations that determine the normal components of v algebraically in terms of the tangential 
components. 

We will show that, under the hypotheses of Theorem [21 the embedding y 0 can be chosen 
so that the tangential subsystem becomes strongly symmetric positive after a fairly simple, 
but carefully chosen, change of variables. Consequently, it follows from the argument given 
in the proof of Theorem |T] that the tangential components of v satisfy the smooth tame 
estimates required to implement a Nash-Moser iteration scheme for the isometric embedding 
system ( 11 . 11 ) . and then the remaining algebraic equations will imply the necessary estimates 
for the normal components of v. Theorem [2] then follows directly from the Nash-Moser 
implicit function theorem (cf. Theorem IA.1D . 

Notation 5.1. We will use the Einstein summation convention for the remainder of this 
paper. 

6. The linearized isometric embedding system and Nash-Moser iteration 

Let fl c K n be a neighborhood of x = 0. Let y 0 : 12 —> R N be a smooth embedding, 
and let g = g l3 dx l dx ] be the metric on 12 induced by the restriction of the Euclidean metric 
on R N to y 0 (12). Linearization of the isometric embedding system (II.ip at the function y 0 
yields the linear PDE system 

( 6 . 1 ) diy 0 ■ djv + dj-yo • d t v = hij , 1 ^ *, j ^ n, 

for the function v : 12 —* R N , where h %3 = gij — gij. 

As described in flj, the linearized system (16.11) can be reformulated as a system of n linear 
PDEs for the n tangential components of v, together with a system of (N — n ) algebraic 
equations for the normal components. To this end, note that, since yo is an embedding, for 
each x e 12 the tangent vectors (diy 0 (x),..., <9 n y 0 (x)} arc linearly independent and span an 
n-dimensional subspace T x c R N . We can therefore decompose the second derivatives of y 0 
as follows: 

(6-2) djy 0 = r-)4y 0 + H ijt 

where, for each 1 ^ i,j ^ n, the vector-valued function H i3 = H 3l : 12 —> R N satisfies 
Hij ■ dkyo = 0 for 1 ^ k ^ n. The functions r*- : 12 —>• K are the Christoffel symbols 
of the metric g tJ , and the quadratic form Hijdx l dx J is the second fundamental form of the 
embedding y 0 . 
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Let S n denote the \n{n + l)-dimensional space of quadratic forms on M n , represented by 
symmetric n x n matrices [sp]. For each x e the vectors 7/^(x) determine a linear map 
i7 x : R N —> S n , given by 

tfx(v) = [<^ii(x),V>]. 

We denote the image by fl x = H X (R. N ). Since the kernel of the map i7 x contains T x , we 
have 


(6.3) 


dimll x < -n(n + 1 ) — n = -n(n 


1 )- 


Definition 6.1. The embedding y 0 : —> R N is called nondegenerate if dimll x = |n(n — 1) 

for all xefl. 


Now, let S* denote the dual space to S n , represented by symmetric matrices [s* J ], with 
the pairing <S* x S n —> M defined for A s S*, H e S n by 

n 

(6.4) {A,H)= Yj A>iH ‘i- 

i,j =1 


Definition 6.2. The annihilator 11^ of the subspace II X cr S n is the subspace of S* dehned 
by 

Il£ = {A e S* : (A, H) = 0 for all H e II X }. 

It follows from equation (I6.3|) that dim 11^ ^ n, with equality for all x e Q if and only if 
yo is nondegenerate. 


Assumption 6.3. Henceforth, we will assume that yo is nondegenerate, and that conse¬ 
quently dim II X = \n{n — 1 ) and dim 11 ^ = n for all xeSl. 

Now, the system (16.ip can be rewritten as follows: 

(6.5) di(djy 0 • v) + dj(diy 0 ■ v) - 2v ■ djy 0 = hq, 1 s$ i, j ^ n. 

Define functions h* and Vjh; by 

Vi = v • diy 0 , 1 ^i<n, 

\7jVi = djVi - T^v k , 1 ^ i, j ^ n- 
then the system (16.51) can be written as 

(6.6) ViVj + WjVi - 2v • H i:j = h i:j , 1 i, j ^ n. 

Since dim 11^ = n, there must exist smooth maps A 1 ,...,A n : Q —» S* such that, for 
each xeH, the matrices A x (x),..., A n (x) comprise a basis of II^. By writing A k = [AW] 
and pairing each of these matrices with equations (16.61) as in (16.41) . we obtain the following 
system of n first order PDE’s for the functions hi, ..., v n : 

(6.7) A kij (ViVj + VjVi) = A kij h ij: k = l,...,n. 

Because A k e S* , we have A k,J = A fcj *, but the component functions A kl do not necessarily 
possess any other symmetries. 

Proposition 6.4. Any solution (hi,..., v n ) : fl —>■ M n to (16.71) uniquely determines a solution 
v : D —» R n to equation (16. lj) .• moreover, v can be determined algebraically from (hi,... ,h n ). 
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Proof. Suppose that v±,... ,v n satisfy (16.71) . and define 

( 6 . 8 ) rjij = hij - WiVj - VjVi, 1 <i,j^n. 


Equation (16.71) implies that [ 77 ^(x)] e II X for each x e Q. Assumption 16.31 implies that // x 
has maximal rank, so that, for each xeO, there exists a unique v(x) e WL N such that 


v(x) ■ diyo(x) = Ui(x), 1 < i < n, 

-2v(x) • H^j (x) = rjjj (x). 1 K:i,j<n. 

Therefore, the map v : ff —> satisfies (16.61) . which in turn is equivalent to (16.11) . □ 


It follows from Proposition 16.41 that, in order to solve the linearized equations (16.11) . it 
suffices to solve equations (16.71) . This system can be written as 



( 6 . 10 ) 


Since A klJ = A kjl and Tf- 


( 6 . 11 ) 


We can write this system in matrix form as follows: For i = 1 let A 1 denote the 

matrix 


- t ' j^lin 


A i = [A kij ] 


nil ' ' ' Anin 


Then the system (16.111) can be written as 


( 6 . 12 ) 


A l diV + Bv = h, 


where 


(6.13) 



V 


Our proof of Theorem [2] is based on the following key result. 

Proposition 6.5. Suppose that the system (I6.12j) is strongly symmetric positive at x = 0. 
Then there exist a neighborhood, fl 0 c ^ of x = 0, an integer (3, and e > 0 such that, for 
any C x metric g on Oo with \\g — g\\p < e, there exists a C x solution y : O 0 —>• ® N to the 
isometric embedding system (13. 

Moreover, the conclusion holds if the system (16.121) becomes strongly symmetric positive 
after performing a change of variables of the form 


(6.14) 


x = 0(x) 


w = 


where (j >: O —► M n is a local diffeomorphism of Pi with </>(0) = 0, and S : fl —* R nxn is a C°° 7 
7i x n matrix-valued function on 0 with S'(O) invertible. 

Proof. First, suppose that the system (I6.12p is strongly symmetric positive at x = 0. The 
argument from the proof of Theorem |T] shows that, under the hypotheses of the proposition, 


there exists a neighborhood f2 0 c= of x = 0 on which the system (I6.12p corresponding to 


the linearization of (II.ip at any function y : flo “^ ~^L N sufficiently close to yo has a solution 
v that satisfies estimates of the form 
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k ^ 0 


for some constants C' k . Then it follows from equation (I 6 . 8 [) that 

\\v\\k < Cl ^||h|| fc+1 + ||h|| 0 ||y - y 0 ||fc +4 +[§]) , 

for some constants C'^. These estimates, together with equations (16. 9 p and Assumption 16.31 
imply (possibly after shrinking f 2 0 ) that 

(6-16) ||v|| fc ^ C k f||h||fc+i + |h|| 0 ||y - y 0 |U +4 +[f]) , k ^ 0 

for some constants Ck- The existence of a solution y : f2 0 ~ N to the system (11.11) then 
follows from Theorem IA.11 just as in the proof of Theorem [TJ 

For the second statement, assume that Qo has been chosen so that the restriction of (j) to f2o 
is smoothly invertible and the matrix S'(x) is invertible for all x e with the determinant 
of 5(x) bounded away from 0. Then it suffices to observe that a change of coordinates of 
the form (16.141) induces linear maps tjjk : H k (Q 0 ) —* H k (cj)(Q 0 )) dehned by 

V’fc(v) = w, 

and that these maps are continuous with continuous inverse. Thus estimates of the form 
(16.151) for the function w imply similar estimates for v, which in turn imply the estimates 
flPSD for v. ^ □ 

Thus it remains to show that, under the hypotheses of Theorem [21 the approximate em¬ 
bedding y 0 : 0 —> ~R n can be chosen so that the linearized system (I6.12p becomes strongly 
symmetric positive at x = 0 after a change of variables of the form (16.141) . 

7. Symmetrization 

The matrices A 1 in the system (I 6 . 12 p are not necessarily symmetric, because the functions 
A k u and A> lk are not necessarily equal. The system (I6.12p can be re-expressed as a symmetric 
system if and only if there exists an invertible n x n matrix C such that the matrices 

CA\ ..., CA n 

are all symmetric, in which case multiplying the system (I6.12p by C results in a symmetric 
system. 

Observe that multiplying (I6.12p by an invertible matrix C is equivalent to replacing the 
basis A 1 ,..., A n for the annihilator 11^ at each point with the alternate basis 

A' 1 = C[A k . 

Moreover, a given basis A 1 ,.. ., A n for 11 ^ will lead to symmetric matrices A 1 , .. . , A n if and 

only if 

(7.1) A kij = A jik - 

i.e., if and only if the coefficients A^ k are symmetric in all their indices. Therefore, in order 
to determine whether the system (16.121) is symmetrizable, it suffices to determine whether 
there exists a basis A k = for 11^ for which the coefficients A kl i are symmetric in all 

their indices. If we choose such a basis for 11^, then we will have 

A k = A k , 

and there will be no need to distinguish between the two. 

Proposition 7.1. When n = 2 or n = 3, the linearized system (I6.12p is symmetrizable. 

20 































Proof. When n = 2, we have N = 3. Choose any smoothly varying basis element ih 3 (x) e II X . 
Consider the 4-dimensional space of all symmetric cubic forms 

A = A ijk -JL-') e S 3 (TR 2 ), 

\dx l dx 3 dx k J 

and for k = 1,2, let A k denote the matrix A k = The annihilator equations 

(A k ,H 3 ) = 0, k = 1,2 

form a system of 2 homogeneous linear equations for the 4 functions A kl] . Thus there must 
be at least a 2 -dimensional solution space at each point x e U, and choosing R(x) to be 
any smoothly varying, nonvanishing element of this space produces a symmetric linearized 
system (16.121) . 

When n = 3, we have N = 6. Choose any smoothly varying basis (LT 4 (x), IL 5 (x), IL 6 (x)) 
for the space ll x . Consider the 10-dimensional space of all symmetric cubic forms 

A = A'* (f- of-of-A S \TV?), 

\OX 1 OX 0 OX K ) 

and for k = 1,2,3, let A k denote the matrix A k = The annihilator equations 

(A k ,H a ) = 0, A: = 1,2,3, a = 4, 5, 6 

form a system of 9 homogeneous linear equations for the 10 functions A kl P Thus there must 
be at least a 1 -dimensional solution space at each point x e U, and choosing R(x) to be 
any smoothly varying, nonvanishing element of this space produces a symmetric linearized 
system (I 6 . 12 K . □ 

Remark 7.2. The result of Proposition 17.11 does not hold for a generic choice of Il x when 
n ^ 4; this is the primary obstruction to applying our methods to the isometric embedding 
problem in higher dimensions. 

For the remainder of this paper, we will restrict to the cases n = 2 and n = 3. We will 
assume that the functions A k ^ are symmetric in all their indices, so that the matrices A 1 in 
the linear system (16.12p are symmetric and may be identified with the matrices A\ We will 
use the convention that Roman indices ( i,j,k , etc.) range from 1 to n, while Greek indices 
(a, ft, 7 , etc.) range from (n + 1) to IV = \n(n + 1). 

8. Compatibility equations and normal forms 

In this section, we will show how the Gauss and Codazzi equations (also called the “com¬ 
patibility equations”) for the embedding yo : 12 —> R N introduce constraints on the values 
of the matrices A 1 (now assumed to be symmetric) and their first derivatives at x = 0, and 
we will show how the matrices A 1 can be put into a simple normal form at the point x = 0. 

Let x = (x 1 ,... , x n ) be local coordinates on centered at x = 0. We will assume that x 
is a normal coordinate system at 0 with respect to the metric g on O, i.e., that Th( 0 ) = 0 
for 1 ^ i,j,k ^ n. We will not , however, assume that g t j(0) = 6ij, because our argument 
will involve a nontrivial GL(n, R) action on the tangent space T 0 M. The specific values of 
gij (0) will not affect our argument, in any case. 

Let g be a real analytic metric on that agrees with g up to order at least /3 (where ft 
is as in Proposition 16.511 at x = 0, and note that this implies that the Riemann curvature 
tensors of g and g agree up to order at least (/3 — 2) at x = 0. By the Cartan-Janet isometric 
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embedding theorem [13], there exists a real analytic isometric embedding (possibly on a 
smaller neighborhood) yo : G —> R w of (Q,g) into R^. 

Let (e n+ i,.. ., ejv) be a smoothly varying orthonormal basis for the normal bundle of the 
embedded submanifold y 0 (G) c R' v , chosen so that 

( 8 . 1 ) Vie a ( 0 ) = 0 

for n + 1 ^ a ^ N and all w e T 0 M , where V 1 denotes the connection on the normal 
bundle induced by the Euclidean connection on . (This condition is the analog for the 
normal bundle of the normal coordinates condition T^-(O) = 0.) Then we can write the 
second fundamental form of yo as 

(8.2) Hijdx 1 o dx j = e a (x) H^dx 1 o dx j 

for scalar-valued functions —» R. 

The embedding y 0 : G — * R w must satisfy the following conditions at x = 0: 

• Metric conditions: 


(8.3) 

(8.4) 


(^iyo • d,yo)|x=o = 9ij( 0), 1 < i, j < n, 

Vi(djy 0 )| x= o = (rjd fc y 0 )| x=0 = 0, 1 < i,j , k ^ n; 

Gauss equations and their first derivatives: 

N 

= Riju( 0 ), 1 ^i,j,k,£^n, 


V 1 ( T-T (y TJ a TJ a U a \ 

1=71+1 
N 

< l L wys - H « H A 


ia=n+l 


x=0 


x=0 


= (d m Rijkt)( 0 ) , 1 ^i,j, k, £, m ^ n, 


where R %3 u denotes the components of the Riemann curvature tensor of ( M,g ). 
• Codazzi equations: 


(8.5) = (<5jffS)|„ 0 = (4ffg)|„ 0 , 1 

This form of the Codazzi equations at x = 0 relics on the normal coordinates condi¬ 
tion T} fc ( 0 ) = 0 and the condition (18.11) on the covariant derivatives of e a . 

Conversely, the Cartan-Janet theorem guarantees that for any choice of real numbers iL“ (0) 
and d k H°j(0) satisfying equations (I8.3p - (I8.5I) . there exists a real analytic isometric embedding 
y 0 : G —> R w of (Q,g) (possibly after shrinking Q) whose second fundamental form agrees 
with the given values up to first order at x = 0 . 


Notation 8.1. Henceforth, we will only be concerned with the values of A k>j . and their 
first derivatives at x = 0. Thus we will use the following notations: 

• will denote the real number //“(0), and H a will denote the matrix [HA]. 

• hfaj will denote the real number dkHfj{ 0), and h% will denote the matrix [h+j]. Note 
that the Codazzi equations (I8.5[) are equivalent to the condition that the are 
fully symmetric in their lower indices. 

• A kl i will denote the real number H fcy (0), and A k will denote the matrix [CT A?? ]. 

• a k + will denote the real number d(,A kl \ 0 ), and a k will denote the matrix [a^]. Note 
that the a k + are fully symmetric in their upper indices, but there are no symmetries 
involving the lower index. 
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• Rijki will denote the real number Rijkt{ 0). Note that the Riju must satisfy the 
symmetries of the Riemann curvature tensor: 

Rijki Rjikt Rijik Rkiij- 

When n = 2, the only nonzero component of R is the Gauss curvature K = 2 ; 

when n = 3, R has 6 nonzero components, represented by R 1212 , -R 2323 , -R 3131 , -R 1223 , 

-^2331, -^3112- 

• fijke,m will denote the real number d m Rijke( 0); when n = 2, we will denote ri 2 i 2 ,i 

and t 1212,2 by k\ and /C 2 , respectively. Note that the must satisfy the same 

symmetries as the Riju in their first four indices, together with the second Bianchi 
identities. When n = 2, the second Bianchi identities are trivial; when n = 3, they 
are represented by the three equations 

( 8 . 6 ) r 2 323,l + f 2331,2 + r 1223,3 = f2331,1 + r 3131,2 + 73112,3 = r 1223,l + ^3112,2 + ^1212,3 = 0. 

The values of Hf-. h“ fc , A kvj . a kl] are constrained by the following relations and are other¬ 
wise arbitrary (apart from the nondegeneracy condition on the HA): 

• Gauss equations: 

N 

(8.7) D-Rnu, 1 < i,j, k, l < n; 

a=n +1 

• Codazzi equations: 

(8.8) h“ jk = h* kj = h“ fc , 1 ^i,j,k^n, n + 1 ^ a < N; 

• Annihilator equations: 

(8.9) A kij H« = 0, 1 < k < n, n + 1 < a < iV; 

• Derivatives of the Gauss equations: 

N 

(8.10) Y, ( H ?k h °im + H ji h ikm ~ H it h< jkm ~ H fk h tim) = r ijU,m , 1 < b 3, k, l, 771 b 71J 

a=n+l 

• Derivatives of the annihilator equations: 

(8.11) A kij h^ j( , + H«af = 0 , n+l^a^N. 

It will be helpful to reduce to the case where the values and A klJ take on relatively 
simple normal forms. To this end, consider a linear transformation of the independent 
variables of the form 


( 8 . 12 ) 


g • x 


with g e GL(n,M.). This transformation induces an analogous action by g on the tangent 
and cotangent spaces T 0 M n and Tq M n , and hence on the tensors 

R = Rijke{dx l a dx 3 ) o ( dx k a dx e ), 

H = e a (x) Hfjdx 1 o dx 3 , 


A 


j\ijk 


8 8 8 \ 

-r O -: O - 

dx l dxi dx k J 


and their covariant derivatives, while preserving the normal coordinates condition T^(O) = 0. 
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$.1. Normal form for n = 2. When n = 2, the subspace IIq cr S-> is spanned by the matrix 


H 3 = 


TT 3 ttZ 

H n H l2 

H 3 H 3 
n 12 n 22 


The nondegeneracy of the embedding y 0 : implies that the matrix H 3 is nonzero. 

Then, by an action of the form (I8.12p . we can arrange that 

K 0" 


(8.13) 


H 3 = 


0 1 


where K is the Gauss curvature of (M, g) at x = 0. The annihilator equations (18. 9 p then 
imply that we can choose 


(8.14) 


A 1 


0 1 
1 0 


A 2 


1 0 
0 —K 


8.2. Normal form for n = 3. When n = 3, II 0 cr dq is the subspace 

II 0 = span (H\H 5 ,H 6 ), 

where, for a = 4, 5, 6 , 


H a = 


Each symmetric matrix H a may also be regarded as representing the quadratic form H^dx 1 dx^ 
e S 2 {T* R 3 ), or equivalently, the quadratic polynomial H^X'X 3 . 

Following [lj, we say that Ho is general if there exists a nonsingular cubic polynomial 
Y = Y ijk X i X 3 X k such that 


H?[ 

TTOL 
11 12 

TTOL 

“31 

TT OL 
** 12 

rra 

n 22 

rra 

“23 

TTOL 

■“31 

TT Ct 

^22 

rra 

“33 


TT ( dr dr dr \ 

n„-sp a n^—j. 


In particular, V must depend on all three variables X 1 , X 2 , X 3 . 
The following classical lemma may be found, e.g., in 


Lemma 8.2. IfYe ^(TqM 3 ) is a nonsingular, homogeneous cubic polynomial, then there 
exists a unique real number a Y —\ and a basis (A" 1 , A" 2 ,X 3 ) o/TqR 3 such that 

Y = (A 1 ) 3 + (A 2 ) 3 + (A 3 ) 3 + 6 cjA 1 A 2 A 3 . 


It follows that, if IIq is general, then by an action of the form (18.121) . we can arrange that 




( 

A 

0 

0 “ 


"0 

0 

a 


"0 

a 

0 “ 

\ 

(8.15) 

II 0 = span 


0 

0 

cr 

5 

0 

1 

0 

5 

a 

0 

0 





0 

a 

0 


a 

0 

0 


0 

0 

1 

/ 


The annihilator equations (18. 9 p then imply that we can choose 


‘- 2 a 

0 

o' 


“o 

0 

r 


"0 

1 

0 

0 

0 

1 

, A 2 = 

0 

- 2 a 

0 

, A 3 = 

1 

0 

0 

0 

1 

0 


1 

0 

0 


0 

0 

- 2 a 
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Meanwhile, the Riemann curvature tensor R may be regarded as a quadratic form on the 
space A 2 (T 0 1R 3 ); as such it is represented by the symmetric matrix 


(8.17) 


R - 


R2323 R2331 Rl223 

R2331 ^3131 R3112 

Rl223 R‘3 112 -^1212 


The only invariant of R under the action (18. 12ft is its signature ( p,q ). The following propo¬ 
sition is a direct consequence of Theorem F in [1]; we will give an independent proof below. 


Proposition 8.3. If R is nonzero, then the Gauss equations (18. 7p have a solution ( H 4 , 
H 5 , H 6 ) whose span II 0 is equivalent under the action (I8.12p to the normal form (j8 . 15jl for 

some a with 0 < |cr| < A. In fact, a may be chosen arbitrarily within this range, the only 

restrictions being that: 

• If the signature of R is (1, 0), then we must have a < 0; 

• If the signature of R is (0,1), then we must have a > 0. 


Proof. Let H 4 , H 5 , H 6 denote the basis 



"l 

0 

o' 


“0 

0 

a 


“0 

O 

o' 

II 

0 

0 

cr 

, H* = 

0 

1 

0 

, H* = 

cr 

0 

0 


0 

0 

0 


a 

0 

0 


0 

0 

1 


for 1I 0 . Then, for a = 4,5,6, let 
for some invertible matrix [7^]. 


H a 

Now, for (3 


7/3 


TpHl> 

■ 4, 5, 6, let 7^ denote the vector 

^ • 


Then it follows from the Gauss equations (18.7p that the corresponding matrix R is given by 


(8.18) R 


(75-Te) - ^ 2 (74-74) 
o- 2 (74 -75) - ^(76'76) 
^(76-74) - ^(75-75) 


^ 2 ( 74 - 75 ) - 0 ^ 76 - 76 ) 
(76-74) - cr 2 (75'75) 

^ 2 (75'7e) - ^(74-74) 


^(76-74) - <7(75-75) 

cr 2 (75'7e) - <7(74-74) 
(74-75) - <7 2 (7 6 -7 6 ) 


It suffices to show by example that, with the sign restrictions given above, the vectors 
74, 75, 76 may be chosen so as to obtain a matrix R of arbitrary nonzero signature. We may 
achieve this as follows: Let 74,75,76 be linearly independent unit vectors in M 3 , oriented 
so that the angle between any pair of these vectors is equal to the same real number 9. 
Geometric constraints require that 0 < 6 < and hence —4 < cos 6 < 1. Denote cos 9 by 
0; then from (I8.18p . we have 


(8.19) 


R 


0 — a 2 cr 2 (j) — a a 2 (j) — a 
a 2 (J) — a 4> — a 2 a 2 f> — a 
a 2 (j) — a cr 2 (j) — a <f> — a 2 


The eigenvalues of the matrix (18. 19p are 

A = 0(1 + 2cr 2 ) — cr(<7 + 2), (1 — cr)(0 + a + cr0), (1 — cx)(0 + a + <70). 
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Therefore, for 0 < a < |, we have 


sgn (R) 


and for — | < a < 0, we have 

sgn(R) 


r ( 0 , 3 ), 

2 < ^ < CT+1> 

(0,1), 

b + 

I b 

II 

(2,1), 

0 / A / cr(tr+2) 

a+l <( P < 

(2,0), 

1 o-((7+2) 

V - l+2cr 2 1 
cr G+ 2 ) < A, < I 

1+2<t 2 ^ 

1(3,0), 


f(0, 

3), 

-\<f< 

cr((J+ 2) 
l + 2o- 2 ’ 

(0, 

2), 

1 _ a(a+ 2) 
V - l+2cr 2 


(1 , 

2), 

O-M 2 ) ^ A 
T+27 3 " < 

b + 

1 b 

V 

(1, 

0), 

b + 

I b 

II 


1(3, 

0), 

~©- 

V 

b| + 

1 b 

< 1. 


A slight perturbation of the vectors 74 , 75,76 will replace the double eigenvalue of R with 
distinct eigenvalues, which will lead to R attaining the remaining possible signatures ((0, 2), 
(1,2), and (1,1) for a > 0 and (1,1), (2,1), and (2, 0) for a < 0) as f> varies. 

□ 


Remark 8.4. It is possible to show that, when R = 0, all nondegenerate solutions (if 4 , 
H 5 , H 6 ) to the Gauss equations (j8.7j) are simultaneously diagonalizable under the action 
(I8.12p and are therefore equivalent to the normal form (j8.15f) with a = 0. The cubic form 
A thus becomes reducible, with the result that the rank of equations (18. lip with respect to 
the variables drops from 27 to 21. This drop in rank is the main obstruction to carrying 

out our construction when R = 0. 


9. Strong symmetric positivity for the system (I6.12p 

In this section we will prove the following theorem, thereby completing the proof of The¬ 
orem El 

Theorem 9.1. Suppose that either n = 2 and K A0, or n = 3 and R A 0. Then the 
linearized isometric embedding system (16 .1 21) can be transformed to a strongly symmetric 
positive system in a neighborhood of x = 0 via a change of variables of the form 

(9.1) x l = x l + v = (I + x l Si) w, 

where c* k = c*. ■ e K and Si ,... S n are constant n x n matrices. 

In order to prove Theorem 19.11 we will show that, when n = 2 or n = 3, for any given real 
numbers R. l3 u and satisfying the necessary symmetries with Riju not all equal to zero, 

there exist real numbers Hf-, A kl f af' J that satisfy equations (I8.7IMI8.11I) . together with 
a change of variables of the form (19. ip . that renders the system (I6.12p strongly symmetric 
positive at x = 0. 

At first glance, the strong symmetric positivity condition might appear impossible to 
achieve: From the expressions (I6.13K and the normal coordinates condition T^-(O) = 0, we 
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have B(0) = 0. Therefore, symmetric positivity for the system (16.12p would require that the 
matrix 

n 

<3o = -XX 

2 = 1 

be positive definite, while strong symmetric positivity would require that each of the diagonal 
sub-blocks (Qi)u = 2 a\ (no sum on i) of Q i (cf. equation (12.61) ) be positive definite. Clearly 
these two conditions are mutually exclusive, and the situation appears hopeless. However, 
it turns out that a change of variables provides some unexpected flexibility: 


Lemma 9.2. Under the change of variables (19.lj) . the symmetric linear system (16. 12)) with 
associated quadratic forms Qo and Q i at x = 0 is transformed to a symmetric system 

(9.2) A%w + Bw = h, 
with associated quadratic form Qo at x = 0 given by 

(9.3) Qo = + c ljA 3 \ 

and the (■ i,j)th block of Q i (cf. equation (12.6)) ) at x = 0 given by 

{Qi)ij = diA\ 0) + djA^O) 


(9.4) 


= a‘+a{- (c‘ t + faA k + SfA’ + A’Si + Sj A‘ + A‘Sj 


"j 1 "'i fajk 1 ~ ikJ -* 1 | “j 

Proof. According to the chain rule, up to first order at x = 0 we have 

d d j _ k d d d j _ k d 

dx 1 dx 1 ikX dxi ’ dx i dx i ikX dxi 

Therefore, at x = 0 we have 


d 


(9.5) 


d 


sg 9 - ( 35 - fa'j^i ) « 7 + 


dx 


dx j 

d 


= (/ + a*Se) ( - O ik x k --w ) + (S) - ^ k x k Sj) w. 


d 


J ( 


v dx i lk dxi 

Substitution of (19.5j) and (19.ip into the linear system (16. 12p yields 

A i ((I + x e S £ ) 


( d _ j __ k d _ 

I —” - fa jjjw 


\dx 


) + (B(I + x l S e ) + A 1 {Si - i k x k S 3 )) w = h. 


Multiply on the left by (/ + arSe) T , collect terms and relabel to obtain the system 
where 

A i = (/ + x e S e ) T (A i - c) k x k A j ) (/ + x e S e ), 

B = (I + x e S e ) T ( B{I + a : e S e ) + A\Si - d^Sf)) , 
h = (/ + x*S e ) T h, 


(9.6) 


and di now represents AA. Finally, computation of 

n 

Qo = B + B J Yi (Qi)w = di# + d 3 A\ 

2=1 

and evaluating at x = 0 yields equations (19. 3 p and (19.4p . □ 

In light of Lemma 19.21 our strategy for proving Theorem 19.11 will be as follows: 

27 













(1) By applying the GL(n,M) action (I8.12H . we may assume that A kl i and //"• are as in 
equations (I8.13p ~ fl8.14p when n = 2 and as in equations (I8.15p ~ (l8.16p when n = 3. 

(2) Identify values of a\'?■* for which the system (I6.12p can be transformed to a strongly 
symmetric positive system (19. 2 p via a change of variables of the form (19.11) . 

(3) Identify values of hhf k that satisfy equations (18.lip for some a\ from Step (2). 

(4) Show that all possible values of satisfy equations (18.101) for some hf jk from 

Step (3). 

(5) Conclude that, for any R^u not all equal to zero and any there exist //"■, 

A klj , hfj k , and cif ' 3 that satisfy equations (I8.7p ~ fl8.lip and for which the system (16.121) 
can be transformed to a strongly symmetric positive system (19.21) . 

Proof of Theorem \9.1[ 

First we give the proof for the case n = 2. We begin by identifying the values of a\'T 2 for 
which we can arrange that 

(9.7) Qo = A/ 2 , Qi = p/ 4 , 


for given real numbers A, /1 > 0, where I 2 and I 4 denote the 2x2 and 4x4 identity matrices, 
respectively. 

By applying the GL( 2) action (18.121) . we can assume that 


Set 


A 1 = 

'0 

1“ 

, A 2 = 

"l 0 

, H 3 = 

~K 

0" 


1 

0 


0 

>1 


0 

1 


0 0 


so that the matrices A 1 , A 2 , D form a basis for <S 2 , and write the matrices S\ and S 2 as 


r n u 

„12~| 


rji 

„12~| 

s i 

S 1 

, S 2 = 

S 2 

S 2 

_21 

^22 


_21 

^22 

L s i 

S 1 J 


L S 2 

S 2 J 


Then, after some computation, equation (19.31) can be written as 

(9.8) Qo = — ( a l + a 2 ) + (c} 1 + c^QA 1 + (cj 2 + c 22 )A 2 , 

and the equations (19.4p can be written as 

(4)n = 2aJ - 2(0^ - 4 8 *)A 1 - 2(4 + Z U 2 (s))A 2 + 4(sl 2 + Ksf)D, 

(g g) (<3 i) 2 2 = 2a 2 - 2(4 + Z 2 2 i(s))A 1 - 2(4 - 2s4aL 2 + AI<(sf - sf)D , 

(4)i2 = a\ + af — (cl 2 + C44 + Zi2i(s))A l — (c 22 + 4 + Z 122 (s))A 2 
+ 2(4 + //(sj 1 + s 21 - sf))D, 


where Zijk(s) represents a linear combination of the s l f whose precise form is irrelevant. 
Regardless of the values of af '-’, we can set 


(Q ijj ) ^jiii2 
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and solve equations (19. 9 p for the variables s;[ 2 , s 2 2 , s 2 2 (from the coefficients of D ) and 
s 22 , s 2 2 , c\ 2 , c 2 2 , c 2 j, c 22 (from the coefficients of A 1 , A 2 ). Note that this solution makes use of 
the assumption that K ^ 0. Then we can set 

Qo = A / 2 

and solve equations (19.81) for the variables c] j and c 22 if and only if the matrix 

a 211 + a 212 + A a 222 + a 222 
a 222 + a 222 a 222 + a 222 + A 

is a linear combination of A 1 and A 2 , which in turn is true if and only if 

(9.10) (a} 22 + a 222 + A) + /\(a 222 + a 2 22 + A) = 0. 

Thus, the strong symmetric positivity condition (19.71) can be realized if and only if the a kl3 
satisfy equation (19.101) . 

The next step is to identify the values of h^ k that satisfy equations (18. lip for some a kl3 
satisfying equation (19.101) . Equations (18.111) may be written in matrix form as 

(9.11) (A k , h^) + (// 3 , a 1 )) = 0. 

The condition (19.101) is equivalent to 

<// 3 ,a} + a 2 > = —(K + 1)A; 
therefore, (19.111) implies that we must have 

(A 2 , h\) + (A 2 , hi) = -<// 3 , a{ + al) = (.K + 1)A, 


CL^ + fl 2 + XI 2 — 


3hfi2 — Khl 22 — (K + 1)A. 


or, equivalently, 

(9.12) 

Conversely, for any values of h 3 j k that satisfy the condition (19.121) . there exist values of a 1 ) 13 
that satisfy the condition (19. 10[) . 

Finally, consider equations (I8.10p . which can be written as 

I<h \ 20 + h\ u = ki, 

(9.13) 22 I 

I\h 222 + h 112 — k 2 - 


The values of h 3 ]k may be chosen arbitrarily, subject only to the condition (j9.12[) : therefore, 
any given values of and k 2 may be realized by an appropriate choice of h 3 )k . 

We conclude that, for any K 0 and any k\. k 2 , there exist solutions h 3 - k and a kl3 to 
equations (j8.7[) — (j8.1ip that satisfy the conditions (I9.10P and (19.121) . and hence the linearized 
system (| 6 . 12 [) can be transformed to a strongly symmetric positive system via a change of 
variables of the form (19.11) . This completes the proof for n = 2. 

Now consider the case n = 3. The argument is essentially the same as for n = 2, but the 
linear algebra requires a bit more effort. We begin by identifying the values of a klJ for which 
we can arrange that 

(9.14) Q 0 = A/ 3 , Qi = /i/ 9 , 


for given real numbers A, p > 0, where I 3 and Ig denote the 3x3 and 9x9 identity matrices, 
respectively. 
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By applying the GL( 3) action (18.121) . we can assume that 


“-2a 

0 

0“ 


‘O 

0 

l" 


"0 

1 

0 

0 

0 

1 

, 7L 2 = 

0 

-2a 

0 

, 7L 3 = 

1 

0 

0 

0 

1 

0 


1 

0 

0 


0 

0 

-2a 


with 0 < \a\ < i. Let H 4 , H 5 , H 6 denote the basis 


(9.15) 



"1 

0 

0“ 


“0 

0 

a 


“0 

a 

0“ 

II 

'tq 

0 

0 

a 

II 

in 

0 

1 

0 

05 

II 

a 

0 

0 


0 

a 

0 


a 

0 

0 


0 

0 

1 


for II 0 . Then, for a = 4, 5, 6 , we can write 


(9.16) 


H a = 7 %H P 


for some invertible matrix [ 7 ^]. Set 


"1 0 0“ 


1 

O 

O 

O 

1_ 


1 

0 

0 

0 

1_ 

0 0 0 

, d 2 = 

0 1 0 

, D 3 = 

0 0 0 

0 0 0 


0 0 0 


0 0 1 


so that the matrices A 1 , A 2 , A 3 , D\, D 2 , D 3 form a basis for S 3 , and write the matrices 
Si, S 2 , S 3 as 



'si 1 


c 13 l 


r s n 

*2 

s 12 

*2 

« 13 1 

*2 


f s n 

*3 

S 12 

b 3 

q 131 

b 3 

Si = 

s 21 

s 22 

e 23 

6 1 

S 2 = 

c 21 

b 2 

s 22 

b 2 

c 23 

b 2 

S 3 = 

s 2i 

b 3 

q22 

b 3 

q 23 

b 3 


C 31 

_ 6 1 

q 32 

6 1 

c 33 

6 1 


,,31 
_ 2 

q 32 

b 2 

q 33 

b 2 


q 31 
_ 3 

q 32 

b 3 

q 33 

b 3 


Then, after some computation, equation (|9.3[) can be written as 

(9.17) Q 0 = —(a} + a 2 + a 3 ) + (c 31 + c 2 2 + c\ z )A l + (c } 2 + c 22 + c 23 )7l 2 + (c } 3 + c 23 + c 33 ) 2 l 3 , 
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and the equations (19.4p can be written as 
(9.18) 

(Oi)ii = — 2(c 31 — sf — s 33 )^! 1 — 2 {c\ 2 + Zn 2 (s))A 2 — 2 (c } 3 + Zu 3 {s))A 3 

+ 4 a(sl 2 + sf - 2sf)D 1 + 4(sf + asf - 2<j 2 s 1 1 3 )D 2 + 4(sf + as 31 - 2 a 2 sf)D 3l 

(Ql )22 = 2a2 — 2 (c^2 + Z 2 ‘ 2 .l{s))A} — 2{c 22 — — S 33 )v4 2 — 2(c23 + Z 223 (s))A 3 

+ 4(s 31 + asf — 2a 2 s 23 )Ili + ia{sf + s 33 — 2sf)D 2 + 4(s2 3 + asf — 2a 2 sl 1 )D 3l 
(Qi )33 = 2a 3 — 2 (c 3 3 + Z 331 {s))A} — 2(c 3 3 + Z 332 (s))^4 2 — 2(c 33 — s^ 1 — s 3 2 )A 3 

+ 4(s 21 + crs 3 3 — 2 a 2 sf)Di + 4(s 3 2 + as 23 — 2a 2 s^)D 2 + 4cr(s 3 1 + sf — 2s 3 3 )H 3 , 

{Q i) 12 = a\ + a 2 — (c\ 2 + c 2 u + Z\ 2 i(s))A l 1 — {cf + c 2 2 + Zi 2 2 (s))A 2 

— (c\ 3 + c 2 3 + Z 123 (s))A 3 + 2 (s 31 + cr(s} 2 + si 2 + s 33 — 2 S 2 1 ) — 2cr 2 s 23 ) D 3 
+ 2 (sf + a(sf + s} 1 + s 33 - 2sf) - 2a 2 si 3 ) D 2 

+ 2 ( S { 3 + S 23 + a(sf + sf) 2 a 2 (sf + sf)) D 3 , 

(Ql)23 = a 3 + a 2 ~ ( C 13 + C 12 + -^231 (s))^ 1 — ( c 23 + C 22 + Z 232 (s))A 2 

— (4 + 4 + W4A 3 + 2 ( S 21 + , 31 + a( 4 3 + 4 2 ) - 2 a 2 (s 32 + , 23 )) D r 

+ 2 (s^ 2 + a{sf + Sg 1 + s 3 3 — 2s 22 ) — 2a 2 sf) D 2 
+ 2 (s 3 + cr(s 2 + s 2 + s 3 — ls 2 ) — 2a s 3 ) L> 3 , 

(Q 1)31 = a 3 + a 3 — (cf + c\ 3 + Z 311 (s))A 1 — ( cf 2 + c\ 3 + Z 312 (s))A 2 

— (c 3 3 + c 33 + Z 3 i 3 (s))A 3 + 2 (s 21 + (t(s} 3 + s 22 + sf — 2sf) — 2a 2 sf) D 3 
+ 2 (sj 2 + s 32 + a(sf + sf) — 2a 2 (sf + s 3 3 )) D 2 

+ 2 {sf + (j(sJ 1 + sf + sf - 2sf) - 2cr 2 s 3 2 ) D 3 , 

where Zijk(s) represents a linear combination of the s l £ whose precise form is irrelevant. 
Regardless of the values of aff we can set 

( Qij) 11 hi 

and solve equations (I9.18P for the variables 

1,1 2,3 3,2 2,2 1,3 3,1 3,3 1,2 2,1 3,1 3,2 2,3 1,2 1,3 3,1 2,1 2,3 3,2 

S 1 ) S 1 > S 1 1 s 2 I s 2 1 s 2 I s 3 I s 3 I s 3 ) S 1 1 s 2 i s 2 i s 2 i s 3 i S 3 ) S 1 i s 3 ) s 3 

(from the coefficients of D±, D 2 , D 3 ) and 

3,3 1,1 2,2 1 1 2 2 3 3 2 1 2 3 1 3 1 2 3 

S 1 > s 2 ) s 3 ) c l,2) c l,3) C l,2) c 2,3) C l,3) C 2,3) C l,l5 C 2,2) c 3,3> C 2,2) c 3,3) C l,l5 C 2,3) C l,3) C l,2 

(from the coefficients of A 1 , A 2 , A 3 ). This solution makes use of the fact that 0 < |cr| < 2 , 
and while the explicit solution is rather complicated, it should be fairly clear that such a 
solution exists for |cr| > 0 sufficiently small. Then we can set 

Qo = 

31 






and solve equations (19.17jl for the variables c} x , c| 2 , c | 3 if and only if the matrix 


a\ + a 2 + O 3 + A /3 — 


ailylg^lia I „U3 I y „112£¥ „122 A „123 „113 „123 


+ On + On 


at 


+ do + On 


at 


+ do + a 


133 


a 


112 


a 


113 


+ a 2 22 + a 
+ ai 23 + a 


123 


+ a; 


+ a 223 + A 


a 


123 


+ do 23 + a 


233 


133 


a 


123 


+ do 23 + a 


233 


+ a; 


+ ajj 33 + A 


is a linear combination of A 1 , A 2 , and A 3 , which in turn is true if and only if 


,111 


(9.19) 


+ a 2 + a 313 + 2 cr(a } 23 + a 223 + a 233 ) + A = 0, 

a } 22 + a 222 + a 223 + 2 a(a } 13 + a 323 + a 333 ) + A = 0, 

a 1 33 + a 233 + a 333 + 2 a(a } 32 + a 322 + a 323 ) + A = 0. 


Thus, the strong symmetric positivity condition (19.1411 can be realized if and only if the a) 1 ' 3 
satisfy equations (19.19j) . 

The next step is to identify the values of h c - - k that satisfy equations (18.11ft for some a) 13 
satisfying equation (19.1911 . Equations (18.1111 may be written in matrix form as 

(9.20) {A\ hf) + (H a , a\) = 0. 

The conditions (19.1911 are equivalent to 

(H a , a\ + a 2 + a 3 ) = —A, a = 4, 5, 6 ; 


therefore, 

<ZP, a\ + a\ + a 3 ) = <7 , a\ + a 2 + a 3 ) = -( 7 “ + 7s + 7e ) A > 
Then equation (19.2011 implies that we must have 
(A 1 , 7i“> + {A 2 , hf) + (A 3 , / 13 ) = —(H a , a{ + a 2 + a 3 ) = (74 + 75 + 76 )^> 


a = 4, 5, 6 . 

a = 4, 5, 6 , 


or, equivalently, 


(9.21) 6/i 123 2ct(/i 111 + /i 222 + /1333) — (74 +75 +76 ct — 4, 5, 6. 

Conversely, for any values of hf jk that satisfy the conditions (19.211) . there exist values of a 1 ) 13 
that satisfy the conditions (19.191) . 

Equations (18.101) are considerably more complicated here than in the n = 2 case. Taking 
the second Bianchi equations into account, there are 15 equations for the 15 components 
t with left-hand sides that are linear functions of the 30 components We will 

regard equations (18. 10ft as defining a linear map G from the 30-dimensional space 77 of 
h?j k values to the 15-dimensional space 77 of rijke, m values; what remains to show is that the 
restriction of G to the 27-dimensional affine subspace defined by equations (19.21)1 is surjective 
onto 77. 

First, observe that we can write equations (18. 10p in matrix form as 

6 

2 G a h a = r, 

a=4 


(9.22) 


32 
































where G a denotes the 15 x 10 matrix 


TTa 

“22 

0 

0 

-277“ 

0 

0 

77“ 

0 

0 

0 

TTa 

“33 

0 

0 

0 

-277“ 

0 

0 

77“ 

0 

0 

0 

0 

0 

TTa 

“33 

TTa 

“23 

0 

0 

TTa 

“12 

0 

Tja 

“31 

TTa 

“23 

0 

0 

TTa 

“31 

TTa 

“12 

0 

0 

0 

0 

-77“ 

0 

0 

0 

TTa 

“23 

TTa 

“22 

0 

TTa 

“31 

0 

0 

Tja 

“12 

0 

TTa 

“33 

0 

0 

0 

^2 #23 

0 

0 

Tja 

“22 

0 

0 

m. i 

0 

TTa 

“22 

0 

0 

-277“ 

0 

0 

0 

0 

0 

0 

TTa 

“23 

0 

-77“ 

Tra 

“31 

0 

0 

TTa 

“12 

0 

TTa 

■“31 

0 

0 

0 

TTa 

“12 

Tja 

“23 

0 

0 

TTa 

“22 

0 

0 

0 

0 

0 

TTa 

“31 

— 77“ 

“33 

0 

TTa 

“12 

TTa 

“23 

0 

0 

m i 

0 

TTa 

“33 

0 

0 

-277“ 

0 

0 

0 

0 

TTa 

“22 

0 

0 

TTa 

“33 

0 

0 

— 277^ 3 

0 

0 

0 

0 

0 

0 

TTa 

“31 

0 

TTa 

“22 

Tja 

“12 

TTa 

“23 

0 

0 

_ H a 
“12 

0 

0 

0 

0 

TTa 

“23 

Tja 

“31 

TTa 

“33 

0 

0 

0 

0 

TTa 

“23 

0 

0 

Tja 

“12 

-77“ 

TTa 

“31 


h° denotes the vector 

y.OL I" l-,OL i-jOL J-iOL Ua UOL UCY Ta 1-.OL 1~.Oi Ta ”1 

n — '‘'222 "333 "ll2 '''311 ""223 "'122 "331 "233 "l23j ) 

and r denotes the vector 

? = [ r 1212,l T 3131,1 r 2331,1 f3112,1 f1223,1 ^2323,2 r 1212,2 f 3112,2 
^1223,2 T 2331,2 r 3131,3 r 2323,3 ^1223,3 r 2331,3 r 3112,3] T - 

Thus the map G is represented by the 15 x 30 matrix 

[G 4 G 5 G 6 ] 


acting on the vector 


[h 4 l 


h = 


h 5 


h 6 


Now, let G denote the restriction of G to the 27-dimensional subspace defined by equations 
(19.21I) . By solving equations (19.21)1 for h“ 23 and substituting into equation 09.22H . we can 
represent G as 

6 

G(h) = G a h a + f 0 , 

a=A 

33 












where G“ denotes the 15 x 9 matrix 


(9.23) G a = 


TTCt 

“22 

0 

0 

-2tf? 2 

0 

0 

i7n 

0 

0 

TTCt 

“33 

0 

0 

0 

-277 3 q i 

0 

0 

77“ 

0 

§*7f& 

1*773“ 

1*^1 

Tja 

"33 

tt a 
"23 

0 

0 

TTCt 

“12 

0 


-5^fi 

"5^fi 

rra 

"31 

tt a 
-"12 

0 

0 

0 

0 

\ aH vi 

1*^12 

\° H ?2 

Tja 

" 23 

tt a 
"22 

0 

TTCt 

■“31 

0 

0 

0 

TTCt 

“33 

0 

0 

0 

2// £; 

0 

0 

TTa 

“22 

0 


0 

tt a 

"22 

0 

0 

-277“ 

0 

0 

\° H V2 


|ai7f 2 

Tja 

"23 

0 

-tffi 

TTa 

"31 

0 

0 


~H ^ 1 — ^aH 22 

~\ pH % 2 

0 

0 

tt a 
"12 

rra 

"23 

0 

0 

\ a Hy . 3 

\ a H23 

\ a ^23 

0 

0 

TTa 

"31 

— 77“ 

"33 

0 

— 77“ 
"12 

0 

0 

H & 

0 

TT a 

“33 

0 

0 

-277 3 “ 

0 

0 

0 

TTCt 

“22 

0 

0 

TTa 

"33 

0 

0 

-277“ 

3 a ^23 

3 *-^23 

3*^23 

0 

0 

TTa 

"31 

0 

TTCt 

“22 

TTCt 

“12 

"5^33 

-\°H%3 

-77 “-§a 77 “ 

0 

0 

0 

0 

TTCX. 

“23 

rra 

“31 

I^ 3 Q 1 



0 

TTCt 

“23 

0 

0 

TTCt 

“12 

- 77 “ 


h a denotes the vector 


t~ol — L a La La La La La La La L a "1 

n — '*222 '*333 '*112 '*311 '*223 '*122 '*331 '*233J > 

and fo is the vector obtained by evaluating G on the vector h with 

K 23 = J(7“ + 7s + 76°)A, a = 4, 5, 6, 
and all other equal to 0. Thus, it suffices to show that the 15 x 27 matrix 


[G 4 G 5 G 6 ] 


has rank 15. 

In order to compute the rank of this matrix, observe that equation (19.161) implies that 


[G 4 G 5 G 6 ] = [ 7 |] [G 4 G 5 G 6 ], 


where G a represents the matrix G a with all entries 77“- replaced by 77,“. Thus, the rank of 

the matrix [G 4 G 5 G 6 ] is equal to the rank of the matrix [G 4 G 5 G 6 ]. We can compute 
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this rank explicitly: Substitution of the expressions (19.15f) for if 4 , H 5 , H & into (19.23ft yields 




0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 a 

— | a —| a — |a 0 0 

0 0 0 a 0 

0 0 0 0 0 

0 10 0 0 

0 0 0 -a 0 

0 0 0 0 0 

|a 2 |a 2 |a 2 0 0 

0 0 10 0 

0 0 0 0 0 

|a 2 |a 2 ^a 2 0 0 

0 0 0 0 0 

0 0 0 0 —a 


0 10 0 

0 0 10 
0 0 0 0 

0 0 0 0 

0 0 0 0 

—2a 000 
0 0 0 0 

-10 0 0 

0 a 0 0 

0 0 0 0 

0 0 0 0 

000 -2a 

0 0 0 0 

0 0 a 0 

0 0 0 -1 


G 5 


1 0 0000000 
0 0 0 0 -2a 0 0 0 0 

|a 2 |a 2 |a 2 00 00 00 

0 0 0 a 0 0 0 0 0 

0 0 00 -1 0 -a 00 

0 0 0000001 
0 0 0100000 
0 0 0 0 0 0 a 0 0 

— ~a -fa —|a 0 0 0 0 0 0 

0 0 0 0 0 a 0 0 0 

0 0 0 0 0 0 0 -2a 0 

0 0 1000000 
0 0 000 -a 0-10 

0 0 000000a 

|a 2 |a 2 |a 2 0 0 00 00 
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G 6 = 


0 

1 

0 

0 




0 

0 


0 

0 

0 

0 




1 

0 


I<7 2 

3 U 


0 

0 

0 

0 

0 


-\ a 


± a 2 
3 u 


0 

0 

0 

0 

0 


0 

0 

0 

0 

F 

0 

0 


v 


0 

0 

0 

0 

0 


a 


a 


-2a 0 0 

0 0 0 


-1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 
a 0 
0 0 
0 0 


1 

0 


0 

0 

0 

0 

0 

0 


0 0 -2a 

0 0 0 


0 a 
0 0 


0 

1 


0 0 
0 0 
0 0 


0 

-1 

0 

0 

0 

0 

0 


0 

0 

—a 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 


0 ' 

0 

0 

0 

0 

0 

0 

0 

0 

—a 

0 

0 

a 

0 

0 


Then a direct computation shows that the matrix [G 4 G 5 G 6 ] has rank 15; for example, 
the submatrix consisting of columns (2, 3, 6, 7, 9,10,12,14,17,18,19,20, 22, 23,25) has deter¬ 
minant equal to —§|a 3 (a — l) 3 (a 2 + a + l) 2 A 0. Therefore, G is surjective onto 77, and any- 
given values of rij k e, m may be realized by an appropriate choice of hf jk . 

We conclude that, for any R =£ 0 and any there exist solutions 77“, A kl i , 7i“ fc , and 

a k y to equations (I8.7p - fl8.lip that satisfy the conditions (19.191) and (19.211) . and hence the 
linearized system (16.12p can be transformed to a strongly symmetric positive system via a 
change of variables of the form (19.ip . This completes the proof for n = 3. □ 


Appendix A. Theorems from analysis 

Theorem A.l (Nash-Moser-Schwartz-Sergeraert). Let E 0 , F 0 be real Banach spaces, arid let 
E k (resp. F k ), k e N, be vector subspaces of E 0 (resp. F 0 ), such that E k+ i a E k (resp., 

F k+ i c= F k ). Let each space E k (resp. F k ) be equipped with a Banach norm || • || fc such that the 

00 00 

inclusions E k+ 1 E k (resp. F k+ 1 F k ) are continuous. Let E * = [~\E k and F .= [^|Ffc 

k =o fc=0 

be given the intersection topology. Moreover, suppose that there exists a family of linear 
“smoothing operators” S(t) : E 0 —> E x , defined for t e R + , satisfying 


00 f 

(1) ||u — S'(t)u||i ^ : |u| 

(2) ||S'(7)u||j < 


(A.l) 


ji 


t e M + , i ^ j, u e Ef 
t e M + , i ^ j, u e Ei , 


where M.\j are positive real constants. 
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Let u 0 e E x ; let D 0 cr E 0 be a neighborhood of u 0 , and let D k = D 0 n E k for k ^ 0. Let 
$ : Dq —i► F 0 be a C 2 map, and suppose that there exists an integer a ^ 0 satisfying the 
following assumptions: 

For any k ^ 0, <3>(A c ) c= F k . 

There exists a constant C such that, for any u e D a and v e E a , 

||$(u + v) - $(u)|| Q ^ C"|jv|| Q , 

||$(u + v) - $(u) - $'(u)v|| a ^ C"||v|£. 

There exist constants C k > 0 with the property that, for any u e D a , there exists a 
continuous linear map R( u) : F a —> E 0 such that, for all h e F a , 

<h'(u) R( u) h = h, 

and for all k ^ 0, ue D k+a , and h e F k+a , 

||R(li) h||fc 5s C k (|h|| k +a "I - ||h|| a ||u uq| k+Oi) • 

Then there exists e > 0 such that, for any f e with 

||f - $(u 0 )||a < e, 

there exists u e D <*, such that 

$(u) = f. 

The proof of this theorem can be found in [25j and |26]. 

Theorem A.2 (Stein). Let LI c M n be a bounded Lipschitz domain. Then there exists a 
linear extension operator 

£ : L\n) — L 1 (W l ) 

satisfying: 

(1) (£/)|q = f; i.e., £ is an extension operator. 

(2) The restriction of £ to W k,p (Ll) is a bounded linear operator 

£ : W k ’ p (Ll) — IU fc ’ p (M n ), l^p^oo, 0^fc<oo. 

The proof of this theorem can be found in [29], 
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